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In this paper, we address the logarithmic corrections to the leading power laws that govern ther- 
modynamic quantities as a second-order phase transition point is approached. For phase transitions 
of spin systems on d-dimensional lattices, such corrections appear at some marginal values of the 
order parameter or space dimension. We present new scaling relations for these exponents. We also 
consider a spin system on a scale-free network which exhibits logarithmic corrections due to the 
specific network properties. To this end, we analyze the phase behavior of a model with coupled 
order parameters on a scale-free network and extract leading and logarithmic correction-to-scaling 
exponents that determine its field- and temperature behavior. Although both non-trivial sets of 
exponents emerge from the correlations in the network structure rather than from the spin fluctu- 
ations they fulfil the respective thermodynamic scaling relations. For the scale-free networks the 
logarithmic corrections appear at marginal values of the node degree distribution exponent. In ad- 
dition we calculate scaling functions, which also exhibit nontrivial dependence on intrinsic network 
properties. 

PACS numbers: 64.60.aq, 64.60.F-, 75.10.-b 



I. INTRODUCTION 

Scaling laws are an intrinsic feature of second-order 
phase transitions. In their leading asymptotics they are 
power laws that govern the behavior of the (singular part) 
of the free energy and of its derivatives in the vicinity 
of the phase-transition point [l|- For a magnetic phase 
transition, the free energy exhibits universal scaling in 
terms of its inherent variables, the reduced temperature 
t = \T — T c \/T c and the magnetic field h. Beside the crit- 
ical exponents universality manifests itself in universal 
amplitude ratios and scaling functions. Moreover, a sys- 
tem defined on a d-dimensional Euclidean space (which 
we will call a lattice hereafter) becomes scale-invariant 
at the critical point. Its correlation length diverges at 
the transition point r = 0, h = and the pair correla- 
tion function changes from an exponential to a power-law 
decay. The leading exponents that govern these scaling 
laws are related by scaling relations. These form a cor- 
nerstone of the modern theory of critical phenomena [l[ . 

Of special interest within this theory of critical phe- 
nomena are those situations in which the aforementioned 
power- laws require logarithmic corrections 0, Q ■ For d- 
dimensional systems, the most prominent examples are 
numerous spin models at their upper critical dimension 
d c and the g-state Potts model in d = 2 dimensions 
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and q = q c = 4 pj. For spin models the logarithmic cor- 
rections appear when the mean-field power-laws observed 
for d > d c turn to non-trivial power-law dependencies at 
d < d c . For the Potts model, the marginal value q c sepa- 
rates two different phase transition scenarios: for q > q c 
the transition is of first order, whereas for q < q c it is 
of second order. Another, more subtle example is the 
d — 2 Ising model with non-magnetic impurities (see e.g. 
Q and references therein). Similar to the leading criti- 
cal exponents, their logarithmic correction counterparts 
have been shown to obey also a set of scaling relations, 
as detailed in Ref. Q . 

The situations discussed above concern systems with 
well-defined Euclidean metrics and, as is clearly seen 
from these examples, the notion of space dimensionality 
is crucial in defining the situation, where the logarith- 
mic corrections to scaling appear. In this paper we want 
to attract attention to a different circumstance where 
critical behavior requires logarithmic corrections to scal- 
ing, namely spin models on networks or random graphs 
0. For Euclidean lattices the space dimension implies a 
given coordination number (2d for the d-dimensional hy- 
percube). For the networks we will consider here, these 
coordination numbers (or degrees) are distributed ac- 
cording to a given degree distribution. This amounts to 
a difference of principal between the origin of logarithmic 
corrections on regular lattices and on such networks. 

The interest to study critical phenomena on complex 
networks is motivated by a number of reasons (Toj both 
of academic and practical nature. Some models on com- 
plex networks may describe exotic phenomena (such as 
opinion formation in a social network fill ]) as well as 
traditional physical objects (e.g., integrated nanoparticlc 
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systems with complex geometry [12j). Real-life complex 
networks are often characterized by a scale-free behavior: 
a power law decay of the node degree distribution 



P(k) 



,-A 



(1) 



Here, P(k) is the probability that an arbitrary chosen 
node of a network has a degree (the number of links at- 
tached to this node) k. The exponent A is crucial in deter- 
mining the critical behavior of different models on com- 
plex networks (see Refs.0, and references therein). 
The general situation is as follows: for small A < A s the 
system is always ordered, only an infinite temperature 
field is able to destroy the order. For large A > A c the 
phase transition is described by the usual mean-field criti- 
cal exponents, whereas systems with intermediate values 
A s < A < A c arc generally described by A-depcndent 
critical exponents. It is the marginal value of A = A c at 
which the logarithmic corrections to scaling appear as has 
been established for a number of classical spin models on 
scale- free networks [l3| . The emergence of these correc- 
tions signals the relevance of correlations between node 
degrees due to the presence of high-degree nodes (hubs) . 
Here, one observes a certain similarity with the critical 
behavior on lattices, where the logarithmic corrections 
appear at the upper critical dimension d c at which the 
trivial mean-field exponents turn to the non-trivial ones 
due to the correlations in thermal fluctuations. 

In the present work we pay special attention to the 
analogy between the role of the upper critical dimension 
d c on a regular lattice and the exponent A c on a complex 
network. To this end, we consider the field and tem- 
perature dependencies of thermodynamic quantities that 
characterize the system in the vicinity of the phase tran- 
sition. The specific example we consider is a system with 
two coupled order parameters on a scale-free network. 
This model is widely used to describe ordering phenom- 
ena in systems with two possible types of ordering. Phys- 
ical examples are given by ferromagnetic and antiferro- 
magnctic, ferroelectric and ferromagnetic, structural and 
magnetic ordering [lj| . In sociophysics applications [ll[ , 
one may think about opinion formation where a coupling 
exists between the preferences for a candidate and a party 
in an election. Recently, we have used a Landau-like ap- 
proach and a mean-field analysis to obtain the phase dia- 
gram of this model on a complex scale- free network [l5| . 
In the present paper we extend this analysis to derive the 
full set of critical exponents that govern the scaling laws 
for the thermodynamic quantities in terms of functions 
of h at r = and of functions of r at h = 0. A special fo- 
cus of our paper is the logarithmic-correction-to-scaling 
behavior. We check the validity of existing relations for 
the logarithmic-correction-to-scaling exponents and fur- 
ther derive new scaling relations for exponents of loga- 
rithmic corrections, for which these relations were so far 
unknown. 



II. CRITICAL EXPONENTS AND 
LOGARITHMIC CORRECTIONS TO SCALING 

The behavior of a system near a second-order phase 
transition is described by a number of critical exponents. 
The magnetization M, susceptibility x an d heat capacity 
Ch at zero external field, h = 0, respectively follow the 
power laws [l(| 



M 



(2) 



(3) 



(4) 



Spatial characteristics of the system, namely, the cor- 
relation length and the correlation function, which are 
connected with the linear size and the spatial dimension 
d, scale with their critical exponents v and r\ correspond- 
ingly. The exponents connected to the spatial structure 
of the lattice are not well defined for the network. At the 
phase transition temperature r = the dependencies of 
the thermodynamic characteristics on the external field 
are also described by a number of critical exponents (l6j 



M ~ h 1 '^ 



C h 



(5) 



(6) 



(7) 



The eight critical exponents listed above depend just on a 
few parameters - spatial dimension, spin dimension and 
symmetries of the model. Therefore, from a knowledge 
of just two of the exponents as well as the dimension, 
any other may be determined. Indeed, the remaining 
six exponents are related via the following four scaling 
relations: 



a + 2/3 + 7 
005-1) 

7c 

a c 
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(8) 
(9) 

(10) 
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For the rf-dimensional lattices, the behavior ©-(O is 
valid from the lower to the upper critical dimension. Be- 
yond the upper critical dimension, the exponents become 
those predicted by the mean field approximation. Just at 
the upper critical dimension one may see modifications 
to the dependencies described above: there appear log- 
arithmic corrections Q. In the absence of an external 
field (h = 0) the scaling behavior at the upper critical 
dimension is 
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X~r 7 |lnT| 7 , 



(13) 



C h ~T- a \\n T \ & , (14) 
while at the critical temperature (r = 0) one finds 

M ~ h 1/S \\nhf, (15) 



X-Zi-^Iln/^, 



(16) 



(17) 



These hatted exponents for the logarithmic corrections 
are also connected via scaling relations, and in Ref. Q 
the following formulae, which are analogous to (|SJ) and 
(El), were derived: 



0(5-l) = <W-7, 



d = 2/3-7. 



(18) 



(19) 



As it was outlined in the Introduction, we are inter- 
ested in scaling laws for the magnetic phase transition on 
networks with, generally speaking, undefined Euclidean 
metrics. Therefore, the exponents we will be interested in 
are those given by Eqs. ©-(El) that do not involve the 
space dimension d. The scaling relations for them are 
given by (|5])- (|lip . However, only two corresponding re- 
lations for the hatted exponents, (fTSj) . (fT^j) . are available 
in the literature so far [3j. Therefore, before we proceed 
further, we derive in the next section the scaling relations 
for the exponents % (jTB"]) and a c (fTT)) that characterize 
logarithmic corrections to the field-strength dependency 



III. NEW SCALING RELATIONS FOR 
LOGARITHMIC CORRECTIONS 

In 0| a Lee- Yang analysis was used to derive relations 
between the logarithmic-correction exponents, which are 
analogous to the conventional scaling relations between 
the leading exponents. Here, these considerations are ex- 
tended to deal with logarithmic corrections to the field- 
dependency of the susceptibility The Lee- Yang analysis 
concerns the zeros of the partition function in the plane 
of complex magnetic field. The locus of such zeros termi- 
nates at the so-called Yang-Lee edge tyl, which is tem- 
perature dependent. Following Q, we account for the 
possible existence of logarithmic corrections to the scal- 
ing of the edge near the phase transition, and write 



In-, 



(20) 



The gap exponents A and A are related to the more 
conventional exponents through the relations Q 



A = /3 + 7, A = /3-7. 



(21) 



In Q , the Gibbs free energy is written as a function of r 
and h as 

/>oo 

$(r, h) = 2Re / In (h - h(r, r))g(r, r)dr, (22) 

in which ft,(r, r) is the locus of Lee- Yang zeros in the 
complex h plane and where g{r, r) is their density. Inte- 
grating by parts yields, for the singular part of the free 
energy, 



$(r, h) 



-2Re 



G(r, t) exp (i(j))dr 
h — r exp (iip) 



(23) 



where G(r, r) = J^ V ^ T \ g(s, r)ds is the cumulative distri- 
bution function for the zeros, the locus of which is as- 
sumed to be h(r,r) = rexp(i^) (the Lee- Yang theorem 
gives 4> = t/2). In contrast to [3j, where h was set to zero 
in (|2"3"]) , the external field is now kept as a variable here 
in order to determine its contribution to scaling near the 
critical point. From Q, the integrated density is 



G(r,r)=xiV(4) 



(24) 



The functional form of I(x) is undetermined here, but our 
considerations shall not require such details. Introducing 
this into (1231), one finds 



$(r» =xr YL J> ( — 



where 



-2Re 



I(x)dx 



y exp I 



(25) 



(26) 



The specific heat is given by the second derivatives of the 
free energy with respect to r, and is 



(27) 



C(T,h) = xr YL T F<t>[ 

\TYL 



Now, from (j!3p and (|2"0"]) . one may express the scaling of 
the specific heat in terms of that of the edge: 



C(r, h) 



r YL A "llnrYLl 1 ^^ 



h 

TYL 



— , (28) 



which may in turn be written as 

o T 2 
. — ~7T T 



C(r,h) 



\\nh\ 



(29) 



Now it is a simple matter to let r — > so that tyl — > 0, 
and the undetermined function JF'^ becomes a constant, 
yielding 



C{h) = h 2 ~i-*\\nh\^ + \ 



(30) 
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From the leading behaviour one recovers (fTTj). The cor- 
rection exponents lead to the new scaling relation 



. ( 7 + 2)A „ 
o-c = 7 h 7, 



which, from (j2"Tj) yields 

( 7 + 2)0-7) 



+ 7- 



(31) 



(32) 



/3 + 7 

Eq. (|10[) and its logarithmic counterpart 

7c = 5. (33) 

are far more trivial to derive and follow from a single 
differentiation of (fTS"]) with respect to /i. The latter two 
equations (|32p and (|33p amount the desired scaling rela- 
tions for a c and %. 



IV. THERMODYNAMIC AL FUNCTIONS OF A 
COUPLED ORDER PARAMETER SYSTEM ON 
A SCALE FREE NETWORK 

In the previous section we obtained new scaling rela- 
tions (|3"2"1) . (|3"3"|) for the logarithmic corrections exponents. 
Together with the formulas (|T8|) . (|19[) they form a com- 
plete set of scaling relations for the correction to scaling 
exponents defined in (fT2"]) -(fTT )) . The validity of relations 
(fT5)l . p^|) for spin models on lattices was subject to a 
thorough check in Ref. Q . There, it was shown that the 
relations hold for all models where the corrections are 
known explicitly. In particular, these include short- and 
long-range interacting 0(n) models at the upper critical 
dimension d c — 4, spin glasses, percolation and the Yang- 
Lee edge problem at d c = 6, lattice animals at d c = 8, 
regular and structurally disordered Ising model at d = 2, 
q-state Potts model at d = 2 and q c = 4 (see ) . Now, 
we will proceed further to perform a similar check for the 
case of critical behavior on scale-free networks. 



A. Temperature dependencies 

As a case study, we will consider a rather common sit- 
uation met in phase transition theory, when a system 
exhibits several types of ordering. This manifests itself 
by the appearance of two coupled scalar order parameters 
denoted by x\, xi- In a microscopic description, such a 
system may be realized as a coupling between two Ising 
models, each of them being characterized by its own or- 
der parameter Xi, or as an XY model with a single ion 
anisotropy. Here we consider the Hamiltonian with a cu- 
bic anisotropy term 



H 



N 



where si and Sj are spins on nodes i and j correspond- 
ingly, J and u are the coupling and anisotropy con- 
stants, the index v numbers the components of the two- 
component vector, si ■ Sj = J2t=i s ».iS».j is a scalar prod- 
uct. The notation denotes the summation over 
all pairs of connected nodes of the network. Note, that 
the Hamiltonian (|34|) is the free energy of an n-vector 
anisotropic cubic model in the case n = 2. The lat- 
ter is obtained from the 0(n) invariant free energy by 
adding invariants of the symmetry group B n of the n- 
dimensional hypercube |17| . 

In Ref. (l5| , thermodynamical properties of such a sys- 
tem at h — were analyzed for a complex scale free net- 
work with a power law node degree distribution exponent 
A, as in (JXJ) . As usual for models on scale- free networks 
[ 1 31 ] , the system remains ordered for any finite tempera- 
ture for A < 3, but it possess a second order phase tran- 
sition with finite T c for A > 3. The phase diagram of the 
system is characterized by two different types of order- 
ing, either along the edges or along the diagonals of the 
square in space the of the order parameter x = {x\, 22}. 
In the first phase, only one order parameter component 
is non-zero (x\ 7^ 0, xi = or x\ = 0, xi 7^ 0), whereas 
x\ =1:2/0 in the second phase. The temperature- 
dependencies of the order parameter, the susceptibili- 
ties and the heat capacities were obtained and the ex- 
ponents ©-([U) determined. The marginal value A c = 5 
was shown to separate two different regimes of the phase 
transition: for A > 5 the exponents attain their classical 
mean field values 



= 1/2, 7 = 1, " = 0, 



(35) 



whereas for 3 < A < 5 two out of the three exponents are 
A-dependent: 

/3 = l/(A-3), 7=1, a = (A-5)/(A-3). (36) 

Another prominent feature found at A = 5 for the tem- 
perature dependencies of the order parameter and of the 
heat capacity is the appearance of logarithmic corrections 
to scaling that have the form given by Eqs. (fT2|) - ([T4| . 
The corresponding correction to scaling exponents were 
found to be UM 



j3= -1/2, 7 = 0, d = -l. 



(37) 



(34) 



To complete the analysis of the phase transition in the 
above model and to access the leading and correction- 
to-scaling exponents ((0)-([7]) and (|T5|) ~- (fT7)) . correspond- 
ingly) that govern this transition, it is necessary to ana- 
lyze the field dependencies of the thermodynamical quan- 
tities at r = 0. 



B. General relations 

The starting point of our analysis are the expressions 
for the free energy considered for different A in (l5j 
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within a Landau-type analysis which was further sup- 
ported by the microscopic treatment of the correspond- 
ing spin Hamiltonian (|34|) . In the following, it will be 
more convenient to work within the (T, x) -ensemble and 
to consider the Hclmholtz free energy F(T, x) related to 
the Gibbs free energy $(T, h), Eq. (|22p. via the Legcndrc 
transform 



F(T, x) = ^{T,h)+x-h. 
For A > 5, the free energy reads [lj| 

b,^, A c 



F{r,x) = -{T- 



4^1 



2^2 



4 1 



(38) 



(39) 



where 



-x\ . Apart from the fact that the param- 



eters a, b 7 c in (|39"]) are A-dependent (see [l5| for explicit 
expressions) , the free energy ([55)1 has the form of a usual 
Landau-type free energy of a system with coupled scalar 
order parameters [14| . Therefore, the network structure 
does not change the critical exponents for A > 5. How- 
ever, for A < 5 the leading terms of the free energy are 
modified [151 ]: 



F(T,x) = 



F(T,x) 



^(T-T c )\xf + ^\x 
4 x 



2 , -|^|4 ln - 



1 



(40) 



A = 5, 



(T-T c )|xf 



a 
2 

-x 1 x 2 |a?'| 



I A — 1 



(41) 



3 < A < 5. 



Note that the functional form of the coefficients a, 6, c 
in Eqs. (|39[) -(j4T j) differ, see [HI for detailed formulas. 
However this explicit form is not important for further 
calculations therefore we keep the same notation for the 
coefficients in ([39 ]l -(|4T )l . 

To analyze the field dependencies of the thermody- 
namic quantities we consider that an external magnetic 
field h is pointing along the order parameter component 
x\, h = {h, 0}, and write the system of equations of state 
as: 



( dF(T,£) \ 
V dx\ J : 



f dF{T,$y 

V dx 2 ' t 



(42) 



(43) 



The stable states are determined from the matrix of sec- 
ond derivatives 



In 



2 F{T,x) 



(44) 



For a given state, the stability condition requires the real 
parts of the eigenvalues of the matrix (|44[) to be positive. 



Note that these eigenvalues are the inverse susceptibil- 
ities, longitudinal Xy" 1 an d transverse Xj 1 correspond- 
ingly. To derive the heat capacity, one needs to obtain 
the entropy of the system. Following the definition 



S(T,x) 



/ dF(T,x) 
A df 



(45) 



one finds the entropy as a function of the temperature 
T and order parameter x. Knowing the dependence of 
the order parameter x on the temperature and external 
field x = x(T, h), obtained from the system of equations 
of state (|4"2"| . (|4"3")l . one finds the entropy as a function of 
the temperature and external magnetic field S = S(T, h). 
Now the heat capacity 



C h = T 



( dS(T, h) \ 

v or )) 



(46) 



completes the calculations. Below, we sketch the results 
obtained for the magnetic field dependencies of the order 
parameter, susceptibilities and of the specific heat for 
different values of A at r = 0. 



C. Exponents for the dependencies on the 
magnetic field at r = 

Taking the expressions for the free energy (|3T))) - (|41[) in 
the equation of state (|42|) and (|43|) one finds the stable 
solutions. It turns out (see Appendix) that there are 
always two stable solutions (xi 7^ 0, X2 — 0) and (xi 7^ 0, 
X2 7^ 0). In any case the quantities considered have the 
leading form © - (JT]). For A > 5 one obtains the mean 
field exponents whereas for 3 < A < 5 we arrive at the 
nontrivial A-dcpcndent exponents given in Table HI 



Q 


P 


7 


5 a c 


7c 


A 


A > 5 
3<A<5 A_| 


1/2 

1 

A-3 


1 
1 


3 

A -2 H 


2/3 

A-3 
A-2 


3/2 

A-2 
A-3 



TABLE I: Critical exponents governing temperature and field 
dependencies of thermodynamic quantities for different values 
of A. 

At the marginal value A = 5 the logarithmic correc- 
tions of the form (fTB"]) - (fl7)) are obtained as summarized 
in Table [TT] completed by the gap exponents calculated 
via Eq. (|2Tjt. 

With the data of Tabic [TT1 at hand, it is straightforward 
to check the validity of the scaling relations for the log- 
arithmic correction to scaling exponents (fT8|) , ([T9|) , (|32|) , 
([33| . Moreover, one can see that they constitute a sepa- 
rate family that differs from other logarithmic correction- 
to-scaling exponents. To this end, we give in the Table [TT1 
the value of the logarithmic correction-to-scaling expo- 
nents that arise for d = 2 Potts model at marginal num- 
ber of spin states q = 4 and for the 0(n)-symmetrical 
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Q 


/9 


7 


5 


d c 


fc 


A 


scale-free network, A = 5 


-1 


-1/2 





-1/3 


-1 


-1/3 


-1/2 


q = 4 Potts, d = 2 


-1 


-1/8 


3/4 


-1/15 


-22/15 


-1/15 


-7/8 


d = 4 0(n) model 


4-n 
n + 8 


3 

n+8 


n + 2 
n+8 


1/3 






l-n 
n+8 



TABLE II: Exponents for the logarithmic corrections to scaling laws that appear for several models: spin model on a scale free 
network at marginal value A = 5 (our results); d — 2 Potts model at marginal spin states number q = 4; and 0(n) symmetric 
model at marginal space dimension d = 4 (see 0] and references therein, fisl]). 



model at marginal space dimension d = 4. Table U also 
gives further evidence of the validity of scaling relations 
(f8l)- (fTTj) for the leading scaling exponents in particular 
of those that involve a c and j c - Values of the latter ex- 
ponents for scale-free networks have not been available 
before. 



V. SCALING FUNCTIONS 



capacity. Since their derivation follows in a similar man- 
ner as the above introduced functions (|4"8")) - (fSTj) we do 
not give their explicit expressions here. 

The formulae given above hold for the single scalar or- 
der parameter system, from which we will start our con- 
sideration. Then the system of 0(n) symmetrical vector 
order parameter to — {toi,TO2} and the system of cou- 
pled order parameters will be analyzed, for which one 
may easily generalize Eqs.(|T7 | - ([5T|) . 



Finally, we consider how the underlying structure in 
the form of a complex network affects the validity of 
the scaling hypothesis, and if the latter is satisfied, we 
will find corresponding scaling functions. The hypoth- 
esis states that the singular part of a thermodynamic 
potential of a system in the vicinity of the critical point 
has the form of a generalized homogeneous function [l9| . 
For the Hclmholtz potential this statement can be math- 
ematically written as [2(J 



F(T 1 m)=T 2 - a f ± (m/^), 



(47) 



where the sign ± corresponds to T > T c or T < T c 
respectively. Provided the homogeneity hypothesis for 
the potential (|47| holds, one arrives at the scaling form 
for the other thermodynamic quantities p^ |. In particu- 
lar, below we will make use of three different equivalent 
representations for the equation of state [HI, Hj] . The 
Widom-Griffiths scaling form of the equation of state is 



h = m 5 h±(T/m 1/0 ), 
with the alternative representation 

h = T 0S H ± {m/r^). 



(48) 



(49) 



The scaling form of the magnetization reads (see also 

Mr- 

m = T H±(h/T 05 ), (50) 
and the isothermal susceptibility may be written as 

XT = r-<x±(h/T^). (51) 



Note that taking temperature derivatives of (|47|) one ar- 
rives at the scaling functions for the entropy and heat 



A. Single order parameter 

For A > 5 the Hclmholtz potential F(t, to) for the 
system with a single order parameter (magnetization) m 
may be obtained from Eq. (|39[) substituting X\ — to, x-i = 
0. Then F(t, m) in dimensionless variables is 



F(T,m) = ±—Tm 



1 2 , 1 4 

m , 



(52) 



where the energy is measured in units of Fq, and the 
magnetization in units of toq: 



F Q =a 2 T?/b, m = \/^. 

It is easy to see that F(t, m) scales as 

F(t, to) = r 2 /± (to/t 1/2 ) , (54) 



where 



/ ± (o = ±ic 2 + Jc 4 . 



(55) 



For A = 5 due to the logarithmic corrections in the 
free energy (f4"!J|). the scaling form defined above fails. We 
refrain from giving a scaling function for this case. 

For 3 < A < 5, F(t, to) as given by Eq.([iT]l for the 
single order parameter system reads 



Ffr.m) = ±-rm 2 + -m^ 1 . 
2 4 



(56) 



Now the dimensionful quantities F and to (|53|) become 
A-depcndcnt: 



(aT c ) 



X — 1 

X-3 



(57) 
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Again, one can recast (|56|) singling out the scaling func- 
tion as: 



F(T,m) 



.l/(A-3h 



(58) 



where the scaling function /±(C) acquires a A- 
dependence, 



MO 



2 4 



(59) 



and we have taken into account the A-dcpcndcncc of the 
heat capacity critical exponent a = (A — 5)/(A — 3), see 
Table |H 

In Fig. [1] we show the dependence of F(t, m)/T 2 ~ a 
on m/rP for different values of A above and below the 
critical temperature. Here and below we will plot the 
corresponding scaling function in the region of positive 
values of h and m. 

Now let us consider the scaling functions h±, H±, fi± 
and x± defined by Eqs. (|4"5)) - ([5Tj) . For the case A > 5 one 
finds from (|5"2")) the equation of state 



m ± rm 



0. 



(60) 



Representing (f6"0"|) as defined by (|48[) we arrive at the scal- 
ing function h±(Q that describes the equation of state 
in the Widom-Griffiths form 



MC) = i±C, 



(61) 



and the scaling function, defined by Eq. (|4"9")l readily fol- 
lows 



#±(C) = C 3 ±C 



(62) 



The dependence of h/r 136 on to/t' 3 , given by the scaling 
function H± in Eq. (f6"2")) is shown in Fig. [2J To get the 
magnetization scaling function (i±, given by Eq. (|50|) . 
we first note that [i± express the dependence of m/r^ 
on h/r 05 , which is inverse to the dependence, given by 
the function H± in Eq. (|49[) . Correspondingly, the scal- 
ing function /i± may be easily plotted by exchanging the 
axes in Fig. [5] as shown in Fig. [3] by black bold lines. 
Three analytic solutions of the equation of state (|60)) for 
m give three branches for the function [i±. Above the 
critical temperature only one branch is real and is pre- 
sented by the black solid line in Fig. [3] The black dotted 
line in Fig. [3] displays the scaling function /i_ below the 
critical temperature T < T c , which corresponds to the 
real solution of Eq. . The above described solutions 
are given by the following formulae 



MO T 2 
MO = — 7, — T 



6 <p±{CV 



where 



MO = (108C+12V/81C2 



± 12 



1/3 



(63) 



(64) 



The scaling function for the susceptibility (|51l) may be 
easily obtained from (|63|) using the relation 



and it reads 



x±(0 = ( 



MO 



x±(C) 



± 



d/i±(C) 
dC ' 



(65) 



36 



M0V ± (0 



9C 



V/81C 2 ± 12' 

(66) 

The dependence of xt/t 7 on h/r 135 , described by the 
scaling function ([65]) is plotted in Fig. [3] by light (orange 
online) lines. 

As we have observed above, the scaling hypothesis for 
the Hclmholtz free energy holds for A > 5 and 3 < A < 5. 
To proceed further and to get the scaling functions h±, 
H±, fi± and x± m the region 3 < A < 5 we first derive 
from (|56p an equation of state, which now has the form 



A- 1 



m x 2 ± rm 



h = 0. 



(67) 



Again, expressing h in terms of m and making use of 
Eq. P5)l we get for the scaling function h± that enters 
the equation of state in Widom-Griffiths form, 



MO 



A- 1 



(68) 



Note, that in the region 3 < A < 5 some of the critical 
indices acquire A dependence. In particular, to get (|68[) 
one should take into account that (3 = 1/(A — 3) and 
5 = A - 2, (Table I|. Comparing Eqs. (J68J) and ^ one 
can see, that for 3 < A < 5 the functional dependence 
of h± on £ does not change, and the only difference is 
the A-dependence of the first term in the right-hand side 
of Eq. This is not the case for the function H±. 

Indeed, from (|6"T)) and (|41?]) we get for this function 



H±(0 = 



A- 1 



C A ~ 2 ±C 



(69) 



Now not only the coefficient but also the leading power 
of this function is A-depcndent. 

There is one more observation which follows from the 
comparison of Eqs. (|69|) and (|62|) that express function 
H± for 3 < A < 5 and A > 5, correspondingly. Eq. (|6"2"]) 
allows for an analytic solution, which enables us in par- 
ticular to find an analytic form for the scaling functions 
M±j X± at ^ > 5 (see Eq. (|M|) and Whereas a simi- 

lar analytic treatment is possible for H± at integer values 
of the power (A — 2) (i.e. for A = 4), it is impossible for 
the general non-integer value of 3 < A < 5. Therefore, 
we make use of the graphic representation to show the 
behavior of fi±, x± at different A in Fig. |U 

Comparing the plots in Fig. 2] one observes a partic- 
ular feature in the behavior of the isothermal suscepti- 
bility scaling functions x±- I n the region 4 < A < 5 
the curve x+ is above the corresponding curve X- f° r au 
values of the argument. Qualitatively this resembles the 
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(a) (b) 

FIG. 1: The dependency of F(T,m)/r 2 ^ a on m/r 13 for values of A between 5 and 3. Plot (a) applies to temperatures above 
T >T C while (b) applies to T < T c 0. 



2 1 




-s: 



m/x 



m/r 



FIG. 2: The dependence of h/r ps on m/r' 3 given for A > 5 
by the scaling function H± (|62|> . above (solid line) and below 
(dotted line) the critical temperature. 



FIG. 3: The scaling functions for magnetization and isother- 
mal magnetic susceptibility for A > 5 above (solid lines) and 
below (dotted lines) the critical temperature. Black curves: 
m/r' 3 , light (orange) curves: isothermal susceptibility scaling 
function x± (color online). 



case A > 5, where the usual Landau theory (j52")l holds. 
In particular, such behavior means that the plot for the 
zero-field isothermal susceptibility x( T ) nas the usual 'A- 
shapc' in the vicinity of r = 0, i.e. the left shoulder of the 
curve is lower than the right one. This situation changes 
with a further decrease of A: first at A = 4 both curves x+ 
and X- coincide, and then, for 3 < A < 4 the curve x+ is 
below X-- Again, for the zero-field isothermal magnetic 
susceptibility x( T ) this would mean that its left shoul- 
der (T < T c ) is above its right one (T > T c ): the usual 
'A-shape' turns to a 'mirror- inversed-A-shape'. The lat- 
ter is closely related to the universal amplitude ratios of 
the ma gne tic susceptibility r + /T_ = A — 3, considered 
in Ref. IIU. 



B. Vector order parameter 

Now let us consider a system with an 0{n) symmet- 
ric vector order parameter being in particular interested 
in the n = 2 case, when fh = {mi, 7712}. The free en- 
ergy F(t, fh) for such a system may be obtained from 
Eqs. ([3U]) - (|4"Tj) by excluding the coupling term. So, for 
A > 5, F(T,rh) in the dimensionless variables is 

F(r,m)=±lr\m\- 2 + \\rh\\ (70) 

where the units of measure are given by Eq. (|53|) . For the 
vector order parameter the free energy scaling function 
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h/f 



ps 



(a) 




h/x 



ps 



(b) 
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(c) 

FIG. 4: The scaling functions for magnetization and isother- 
mal magnetic susceptibility for A = 4.8 (a), A = 4.0 (b), 
A = 3.1 (c) above (solid lines) and below (dotted lines) the 
critical temperature. Black curves: m/r 13 , light (orange) 
curves: isothermal susceptibility scaling function \±- 



/±(C) is defined by the equation 

F(T,m)=T 2 - a f ± (rh/Tp), 

and reads 

f±(0 = ±\\C\ 2 + \\(\ 4 - 



(71) 



(72) 



Taking that the magnetic field points along the mi com- 
ponent, one finds that the stable state of the system re- 
quires 



and 



m\ ^ Tm i — h = 0, 



m 2 = 0, 



(73) 



(74) 



as follows from the system of equations of state (|42[) . 
(|4"3"| . Since the field points only along the first com- 
ponent of the magnetization, we get that /i2± = and 
Mi± = M±- The scaling functions h±, H± and fj,±± = f/,± 
coincide with the corresponding scaling functions for the 
single scalar order parameter, and are given by Eqs. (|61[) 
(fFJ3")) where m is to be replaced by mi. Nevertheless, one 
needs to consider two different response functions, that 
describe the reaction of the system on an external mag- 
netic field: the longitudinal and transverse susceptibili- 
ties, which scale as 



X|| 



X-L 



{h/r»). 



(75) 



(76) 



Using the definition for \\\, X± ( sec Eq. (|44")l and below) 
and substituting solutions (f73|) . (|74|) into corresponding 
derivatives of the free energy ([70)) we arrive at the scaling 
functions 



xn±(0= ±i + 3(m±(C)) : 



x±±(C) 



±1 



(m±(C)) : 



(77) 



(78) 



with /i±(() given by Eq. (|fJ3")) . Note that the existence 
of the continuous 0(n) symmetry of the free energy at 
h = allows also to present the transverse susceptibility 
ratio of mi and h at arbitrary h: 



X± 



~~h 



(79) 



Let us repeat the above calculations for 3 < A < 5. 
Excluding the coupling term from Eq. f4"Tj) we get for 
the 0(n) symmetric free energy: 



F(r,m) = ±ir|m| 2 + i|m| 



A-l 



(80) 
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and the scaling function f± follows 

f±(t) = ±lK\ 2 + \\{\ x - 1 - (81) 

The units of measure are given by Eq. (|57|) . Again, the 
system of equations of state (|4*2"j) , requires 



m 2 = 0, 



(82) 



and the equation for mi coincides with Eq. (|67[) for a 
scalar order parameter to: 



A - 1 



m^ 2 ± mi - h = 0. 



(83) 



Properly, the scaling functions h±, H± and fi±± = fj,± 
are equal to the corresponding scaling functions for the 
single scalar order parameter. The susceptibilities scaling 
functions \\\± an d X±± follow 



xn±(0 



±i 



X-L±(0 = 



^«(, 1± (C)) A - 3 



±i + ^(m 1± (C)) a - 3 



(84) 



(85) 



Again, as for the case A > 5 we note that the transverse 
susceptibility \± also for 3 < A < 5 can be recast as the 
Eq. (J79]). 



Mi± = M± coincide with the single scalar order parameter 
given in Eqs. ([oT jl -([S3" jl . For the type of ordering consid- 
ered, [mi, 0], the susceptibility scaling functions read: 



X ||±(C)= [±1 + 3( M± (0) 

xx ± (C) = [±i + ^(m±(0) 



(89) 
(90) 



with n±(() given by Eq. ([6"3")l . 

For the second type of ordering, [m 1 ,m 2 ], the system 
of equations of state ([A.ip . (|A.2[) may be written as 



m\ i ™i — h = 0, 

2 _ 2b+c 2 _ 
2 — 2fc 1 ' 



where 



26 



46 2 



46+ c 



c(46 + c) 



(91) 



(92) 



The coefficient in front of r in Eq. ([9"2")) is positive due 
to the stability condition in the vicinity of the critical 
point at zero magnetic field (h = 0) [ll| as well as at the 
critical temperature (r = 0) with applied field (h ^ 0), 
Eq. (TOl) . 

The Widom- Griffiths form of the equation of state (|4"5)l 
reads 



= mfh±(f / m 1 



1/2, 



(93) 



C. Coupled order parameters 

For the model we consider in this paper, the vector m 
has two components which corresponds to two coupled 
scalar order parameters mi and to 2 . For the case A > 5 
the free energy is given by Eq. ([39j) and the corresponding 
scaling function follows 



/ ± (C)=±5lC? 



7ic1 4 + ^Ci 2 c!, 



(86) 



where the units of measure are given by Eq. (|53l) . 

Let us recall, that the presence of the coupling between 
the order parameters mi and m 2 leads to two possible 
stable states of the system. Besides the state 

[mi,0], with mi 7^ 0,to,2 = 0, (87) 

the system is characterized by an additional stable state: 

[mi, TO2], with Tti\ 7^ 0, m 2 7^ 0. 



with h±((), given by Eq. (|6ip . The equation of state, 
given by (|4T)|) is then 



h = ¥ 3 / 2 H ± (m 1 /f 1 ^), 



(94) 



with H±(() given by Ea.([62"]). 

The order parameters mi and m2, obtained from the 
system (|91[) may be conveniently presented in the scaling 
form 



mi = f 1/2 ^i ± 



m 2 



(Vr 3/2 ), 
(h/f 3/2 )- 



(95) 



The scaling function /ii±(C) coincides with the function 
A*±(C) m (p3"l> of the single order parameter system, and 
for /X2±(C) we find: 



mL(C) 



26 



26 



c 2 Ah 
-Mi±(C) T 



26 



(96) 



In the ordered state [mi, 0], where only mi reflects the 
system magnetization, the scaling functions h±, H± and 



The susceptibilities follow the scaling form (J75J), ([76 
where X||±(C) an d X±±(C) ar e 
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Here we used the notations /ii = /Ui±(C), M2 = A f 2±(C)- 

For the case 3 < A < 5 the Hclmholtz potential (pPI 
follows the scaling form (|7T|) where the scaling function 



i|c?- 1 + ^c 1 2 ciici A - 5 



(99) 



becomes functionally A-dcpendent. 

Similarly as for A > 5, for 3 < A < 5 the system of 
equations of state (|A.30[) , (|A.31[) permits two types of or- 
dering at nonzero external field h ^ 0: [mi, 0] (Eq. (|87p) 
and [mi,m 2 ] (Eq. ([88])). 

For the case [mi, 0] where only mi depends on the ex- 
ternal field, the scaling functions h± and H± coincide 
with the scaling functions given by Eqs. ([59^1 . Simi- 
larly, as for the single order parameter, the analytic treat- 
ments for the magnetization and the isothermal suscepti- 
bilities are impossible for general non-integer A, for which 
3 < A < 5, whereas their shape for different A is shown 
in Fig. [4] Note however, that the isothermal susceptibil- 
ity scaling functions, defined by Eqs. (J7SJ), (|76[) may be 
analytically expressed through the magnetization scaling 
function fii± as 

x,i ± (C) = (±i + (A " 1) i A " 2) [Mi ± (C)] A - 3 )" 1 , (ioo) 



xx±(C) = (±i 



A-l 



c 

26 



[Mi±(0] 



A-3 



(101) 



For the ordering [mi, 7712], the equations of state 
(|A.30p . (|A.31[) do not allow to obtain analytical solu- 
tions for mi and m-2 in the general case for arbitrary 
non-integer A. Nevertheless, these equations enable a 
confirmation of the scaling of the magnetization. Indeed, 
substituting 



mi = r 1 /^- 



■3) 



(ft/- 



_(A-2)/(A-3)> 



(102) 



m 2 = T^-Vfi^h/T^-V/^-V) (103) 



into Eqs. (|A.30[) . (|A.31[) one obtains the following equa- 
tion for the scaling functions [i\± and //2± : 

^i ± (C)(^ ± (0-^i ± (C))|m±(C)I a - 5 =C, (104) 

where the relation between /ii±(C) = Mi± an d M2±(C) = 
H2± is as follows 



± \n± 



A-7 1 ^ 1 1 .-? 1 4 ■ c 



— iM±r + 26^?±Im±I 2 



( A ~ 5 ) C 2 2 n 
17 Ml±^2± - U - 



(105) 



The observed dependence of Eqs. (|104[) , (|105[) on a sin- 
gle variable C = h/r^^ 2 ^^^ serves as evidence of the 
validity of the scaling hypothesis. In a similar way, the 
scaling for h, x\\ an d X± rnay be confirmed. 



VI. CONCLUSIONS 

Although second order phase transitions take place 
only under certain conditions, e.g. r = 0, h = for 
magnetic systems, the singular part of the free energy 
and its thermodynamic derivatives are described by func- 
tions characterized by scaling properties nearby. Such 
properties are found experimentally by measurements at 
various points in temperature and small values of the 
field. Introducing appropriate scaling fields, data col- 
lapse to universal scaling functions is used to determine 
the exponents. On the other hand, logarithmic correc- 
tions present in the temperature and field dependence of 
quantities described otherwise by power laws may disturb 
this collapse. Therefore we have on the one hand cal- 
culated several scaling functions for ranges the network 
parameter A where power laws arc valid and found the 
characteristic dependence of these functions on A even 
when the exponents are independent of A. On the other 
hand we have determined the possible logarithmic cor- 
rections at certain borderline values of A. 

It is remarkable that already within mean field theory 
such logarithmic corrections arise. They are attributed 
to correlations due to the network properties rather than 
fluctuations (which are absent in mean field) of the spin 
properties themselves. However we note that an essential 
assumption in the mean field treatment was the additive 
contribution of the individual nodes of the same degree 
to the network free energy. 

In this paper, we gave a comprehensive description of 
the temperature and field behavior of a system with two 
coupled scalar order parameters on a scale-free network 
in the vicinity of the second order phase transition point. 
Special attention has been paid to the appearance of the 
logarithmic corrections to scaling. For magnetic systems 
on d-dimensional lattices, such behavior arises due to the 
order parameter fluctuations that tend to be strongly cor- 
related in the vicinity of the critical point. It is the space 
dimension d that is definitive for the relevance of such 
fluctuations. The scale-free networks we consider here 
are not characterized by the Euclidean metrics and the 



12 



space dimension. Instead, it is the node degree distri- 
bution function exponent A ([1]) that brings about cor- 
relations present in a network due to its internal struc- 
ture. As A decreases and the node degree distribution 
becomes more and more fat-tailed, the relative number 
of the high-degree nodes (hubs) increases and leads to 
non-trivial critical behavior. This can be related to the 
fact that below A = 5 the fourth moment of the degree 
distribution diverges and below A = 3 the second mo- 
ment ceases to exist. First, for A c = 5, the non-trivial 
dependencies appear. With further decrease of A, for 
A < 3, the systems appears to be ordered at any finite 
temperature. 

One observes a certain formal similarity between the 
behavior of spin systems on d-dimcnsional lattices and 
on scale-free networks with exponent A. Both at d = d c 
and at A = A c the logarithmic corrections to scaling are 
precursors of the change in the critical behavior. This 
formal similarity is further pronounced in a more subtle 
way: the correction exponents found by us for the scale- 
free networks, although numerically different from those 
found on the lattices obey the same scaling rela- 

tions. In the present work we derived two new scaling 
relations for logarithmic correction exponents (|3^|) and 
(f3"3")) . Together with previously found [3| relations (fIB"]). 
(fT9f they form a complete set of scaling relations for log- 
arithmic corrections. 



a(T - T c )x 2 + bx 2 \x\ 2 + -cx\x 2 = 0. (A.2) 

At t = one finds two different solutions, which we now 
detail. 

Solution x\ ^ 0, X2 = 0. In the case of vanishing x 2 
this solution is 



1 = ^^ 3 , X 2 =0. 



(A.3) 



This solution exists and satisfies the stability conditions 
if 

6 > 0, c> -2b. (A.4) 
The longitudinal and transverse susceptibilities follow 

1 



X|1 ~ 3&V3 
2b 2 / 3 



A J 



h~ 2 l 3 . 



2b + c 

The heat capacity at the critical point reads 



Solution xi 0, x 2 ^ 0. For x 2 non-vanishing 



(A.5) 
(A.6) 

(A.7) 
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Appendix A 

In this Appendix we give the deviation of the expo- 
nents governing the field dependencies of thermodynamic 
quantities at r = for different values of A. 



46 



c(46 + c) 



7 V/ 3 , x 2 = ±J-^^x 1 . 



(A.8) 

Note, that the ratio between x\ and x 2 does not depend 
on the strength of the field and depends only on the ratio 
c/b. This solution exists and satisfies stability conditions 
at 



b > 0, -46 < c < -2b 
The susceptibilities read 

"c(46 + c)l 2 / 3 



X± 



8b 

> + c) + VZ 



46 

c(46 + c)v 



ft" 2 / 3 , 



46 



where 



£ = c 2 (86 + c) 2 - 486c(26 + c)(46 + c). 
The heat capacity follows: 



(A.9) 



(A.10) 
h- 2 ' 3 , (A.ll) 



(A.12) 



1. Case A > 5 

In this case the free energy follows (f39|) and the system 
of equations of state reads: 



1 



a(T — T c )xi + bx± |x| 2 + -cx±x 2 = h, 



(A.l) 



C h 



3 (46 + c) 



(A.13) 



From the above solutions we conclusion that the expo- 
nents defined in formulas ©-([7]) arc 



S = 3, 



7c 



a r = 0. 



(A.14) 
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2. Case A = 5 



Given the free energy (|40|) , the system of equations of 
state reads 

a(T - T c )x! + fccilfl 2 In -4- - ^xilxl 2 (A.15) 

|x| 4 

1 3 2 

2i„ 1 "A _ ^ 



^ 1 37*-) In , 
2 2 |2| 4 |x| 2 



o(T - T c )x 2 + bx 2 \x\ 2 In ^ - ^x 2 |x| 2 (A.16) 

x 4 



2„.3 



1 C X^X 



+ 2 XlX2ln M"i"N 2 



0. 



Note that this is a system of transcendent equations and 
one may estimate the solution at weak external field h — > 
0. At t = there exist two solutions of the system (|A.15|) . 

(OS) . 

Solution xi 7^ 0, x 2 = 0. The solution for vanishing 
X2 is 



Xl 



3\l/3 &V3 



6/ (-ln/i) 1 ^ 



x 2 = 0. (A.17) 



This solution exists and satisfies the stability conditions 
if 



6 > 0, c> -26. 
The susceptibilities follow 

xi, = (MyWh-Wi-hih)- 1 ' 3 

2b /6\-i/3 



(A.18) 
(A.19) 



26 + cV3, 
The heat capacity reads 



h- 2/3 (-\nh)~ 1/3 . (A.20) 



Ch = -jT c {- \nh)-\ 



(A.21) 



Solution xi 7^ 0, X2 7^ 0. For nonzero x 2 , the solution 
is of the form 



xi 



2b + c 

This solution exists and satisfies stability conditions if 

b > 0, -46 < c < -26. (A.23) 
The susceptibilities follow 

X|| =xf^ 2/3 (-lnM- 1/3 , (A.24) 



X±=xl J ^ 2/3 (-lnM- 1/3 



where 



(A.25) 



(~2^) 5/3(46 + c)2/3 ' (A - 26) 



X-/ = - ^^ , f-^ 5/3 (46 + c)V3, (A , 7) 



c(86 + c) + V£V 26/ 
and £ is defined by Eq. (|A.12j) . The heat capacity is 
2a 2 



C h 



46 



-T c {-\nhy\ 



(A.28) 



Comparing the obtained solutions with the definition 
of the logarithmic corrections to scaling exponents (TH 
(IT71) . we arrive at 



1 



1 



3, 7c = -3, «Sc = -l. (A.29) 



3. Case 3 < A < 5 

The system with the free energy (|4"Tj) is described by 
the following equations of state 



a(T - T c )xi 



A-l 



6xi|x| A 3 + -cxiX2|x| A_5 



+ ^^cx\x\\x\ x - 7 = h, (A.30) 



a(T-T c )x 2 + — — 6x 2 |x|" 



4 

A- 5 



A-3 1 ± r ~2 — I^|A-5 

2 



cx 2 x^|x| A - 7 = 0. (A.31) 



The system of equations (|A.30J) . (|A.31|) has two solu- 
tions. 

Solution xi 7^ 0, x 2 = 0. For the vanishing x 2 this 
solution is 

This solution exists and satisfies the stability conditions 
if 



6 > 0, c > - - „ 1 b. 



(A.33) 



The longitudinal and transverse susceptibilities read 
X|| =T 1 -J7T A TTz) 1/{X ' 2) h~ {X - 3)/{X - 2 \ (A.34) 



A-2V(A- 1)6 



X± 



(A-l)6 N(A-3)/(A-2) (a _ 3)/(a _ 2) 



(A - 1)6 + 2c V 4 
The heat capacity 



(A.35) 



C,, = 



(7T^) 3/(A " 2) ^ (5 - A)/(A - 2) . (A.36) 



A - 2 V(A 



Solution X! 7^ 0, x 2 7^ 0. For the non-vanishing x 2 
the solution is 
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Xl 



4(1 + H 2 ) 



2\(7-\)/2 



l/(A-2) 



h l / {x - 2 \ x 2 = ±nx u 



where 



The solution exists and is stable if 



(A - 1)6(1 + m 2 ) 2 + 2c A j 2 (1 + fi 2 ) - (5 - A)c/u 2 

2 _ -2(A - 1)6 — (A - 3)c + V(A - 3) 2 c 2 - 4(A - 1)(5 - A)6c 



2(A- 1)6 



6 > 0, -46 < c < - ^ 1 6. 



(A.37) 



(A.38) 



(A.39) 



The susceptibilities follow 



X\\,± 



where 



(A - 1)6(1 + (i 2 ) 2 + 2c M 2 (l + /i ) + (A - 5)c^ 



(A-3)/(A-2) 



{(A - 1) 2 6 + 2c}(l + m 2 ) 3 - (A + 3)(5 - A) CA i 2 (l + M 2 ) ± 



4 1/(A-2) 2(1 + /x 2 ) ^±| y/l -(A-3)/(A-2) ) (A>40) 



D = (1 - ^ 2 ) 2 {[(A - 1)(A - 3)6- 2c](l - /it 2 ) 2 + (5 - A)(7 - A)c/i 2 |" 
+4^ 2 {[(A - 1)6 + 2c] (A - 3)(1 + a^ 2 ) 2 + (5 - A)(7 - A)c/i 2 } 2 

Finally, the heat capacity is 



Ch = 



4(l + /i 2 )( 7 - A )/ 2 



(A - 1)6(1 + ^ 2 ) 2 + 2c/i 2 (l + fi 2 ) + (A - 5) CA t 2 

I 



3 /^- 2 )l±M! a 2 T ,( 5 -A)/(A-2). 



A - 2 



(A.41) 
(A.42) 

(A.43) 



The above results give the following values of the critical 
exponents: 

6 = A- 2, 7c = ^|, «c=xr|- ( A - 44 ) 
Note that in this case all leading exponents that govern 



the field dependencies at r = are A-dependent. 

Summarizing results of Ref. [Hf, Eqs. (|35 |) -(|37 | . and 
those obtained in this section, Eqs. (|A.14|) . (|A.29[) , and 
(|A.44I) we give the values of the leading and correction 
to scaling exponents in Tables HI UU completing them by 
the gap exponents A, A calculated via Eqs. (f2~l"T) . 
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